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three dimensional unsteady lifting surface theory 

IN THE SUBSONIC RANGE 
H. G. Kiissner** 


1* Introduction 


/41.2 * 


The aeroelastic problems of the modern airplane and aircraft 
are so numerous and so complex that a purely experimental 
treatment is bound to fail due to the multiplicity of the 
necessary parameter variations. The closer one comes to the 
boundaries of technology, the more necessary it becomes to apply 
mathematical and theoretical methods to solving a problem, even 
when making simplified assumptions. For this reason unsteady 
i ting surface theories (in spite of their recognized 
deficiencies) are used to a large extent in aeroelastic 
investigations of new airplane designs. In the leading airplane 
works they avoid purely theoretical or experimental directional 

at i°n fl “tter proofs, because every known process still has 
great deficiencies and because the only sufficient protection 
against directional oscillation flutter to date has been 
comparative investigations. 


Since this pressing need exists, much work has been put into 
the development of unsteady lifting surface theory in the last 30 
years. The two-dimensional problem has been solved completely 
for all Mach numbers. Many numerical tables have been calculated 
or the use of these solutions. In the case of the three- /41 3 
dimensional problem, that is. in the case of lifting surfaces of 
mite Wingspan, only minimum progress has been made to date 
because the Kut ta-condi t ion in the subsonic range is a difficult 
boundary condition. First, I should like to present you with a 


♦Numbers in the margin indicate pagination in the original. 

♦♦Max Planck-Inst Itute of Flow Research. Gottingen (West Germany) 
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survey of the general three-dimensional 
assessment of same. Then I should like 
investigations in the subsonic range. 


problem and solution 
to report on my recent 


2 • Basic Ass ump tions of the Lifting Surface Theory 

The question how much stress is placed upon a winged aircraft 
with a sharp trailing edge when in flight could not be answered 
in the 19th century. Classical hydrodynamics supplied the 
paradoxical answer that, under ideal fluidity and stationary 
flow, stress is zero. Only at accelerated motion do stresses 
varying from zero — the so-called Kelvin impulse — occur. 


On the basis of Prandtl’s boundary layer theory, Kutta first 
conceived the idea of considering the viscosity of actual 
fluidity phenomenologically. He did this by assuming the smooth 
flow of the trailing edge, while otherwise figuring with ideal 
fluidity. Kutta’s two-dimensional profile theory offered the 
first solutions which were also physically acceptable. 


For three-dimensional flow, that is, for wings of finite 
wingspan, the Kutta method of conformal representation fails. 

For this reason also, simplified assumptions had to be made. It 
was necessary to limit oneself to infinitely small / 41.4 

disturbances, by which the problem is linearized. Accordingly, 
the thickness of the wing must be infinitely small, so that one 
can even talk of lifting surfaces. 


At the trailing edge of this lifting surface the flow should 
stream away smoothly. The sharp front edge of the lifting 
surface, however, causes the flow to go around it, which leads to 
infinitely large disturbance speeds and pressures. As long as 
this specific pressure is integrable, we can nonetheless 
calculate total reactions upon the lifting surface, that is, upon 
lift, linear and higher torques. 
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For this reason, a lifting surface theory can be a theory 
only approaching the calculation of aerodynamic coefficients. It 
is desired, nonetheless, that under the given circumstances the 
mathematical boundary value problem of the lifting surface theory 
be solved exactly and that no further approximate assumptions be 
made. For the danger exists that accumulating approximate 
assumptions will result in loss of touch with physical reality. 


First, we consider the boundary conditions of lifting surface 
theory. We shall assume frictionless compressible fluidity. The 
positive x-axis points in the direction of airspeed v, the y-axis 
in the direction of the wingspan (See Figure 1). 


Further, t repres 
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Figure 1. Moved lifting surface in the plane z = 0. 


/41.5 
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Further we assume that the lifting surface remains even up to 
infinitely small deviations and lies in the plane z = 0 (See 
Illustration 1.) 

Because we have linearized our problem, the Euler equations 
become 


^ <J>(s,x,y,z) + p( 3 ,x ,y , z ) * 0 , (2) 

~ P 2 (VS ~ (e.x.y.z) = 0 . 

We take the vertical component of velocity at the lifting 
surface, the so-called down current, as given. With the 
exception of the leading edge, it is equal to the gradient of the 
velocity potential perpendicular to the lifting surface 


w(s,x,y) = $ 2 (s,x,y,0) 


(4) 


Down current is determined either by the small 
deformations of the lifting surface or given by 
squall area. According to the Kutta condition, 
pressure of the lifting surface at the trailing 
always be p = 0. With that, our boundary value 
completely determined. 


/41.6 

an atmospheric 
disturbance 
edge should 
example is 


3 . Transformation of the Euler Equations. 


We want to simplify the solution of this task 
carrying out two transformations. First we shall 
Galilei-Lorentz transformation of coordinates and 


by first 
perform a 
t ime : 
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X’ 


(5) 


= x ; y’ = y/Tp 2 ; o’ = s(1-p 2 ) - p 2 x , 

2 - = s /l-p 2 ; f = t(1-p 2 ) - p 2 x . 

Using B = 0.7 as an example, the given lifting surface countour 
reverts to the dotted line in Figure 1, which we shall continue 
to consider. 

Further we perform a Laplace transformation; this changes the 
object function F(s’) which disappears with s’< 0 into the 
diagram function 


£T{ F(s’ )] w*} ■ - J F(s' ) exp ( -u*s ' ) ds' . 
o 

In this equation co* is a complex number. 

By means of these two transformations (5) and (6), the Euler 
equations (2) and (3) become 


(u* - <J>*(x,y,z) + P*(x,y,z) *= 0 

( x 2 -V 2 ) 0*(x,y ,z) = 0 


(7) 

( 8 ) 


The transformation lines have been left out here. The reduced 
frequencies are /41 . 7 


U)» 


1 i) 1 


X 


ll) 1 

c(l-p 5 ) 


= p 




(9) 


In the case of harmonic vibrations, is the purely imaginary 
circuit frequency of the object function. 


7 



4. 


Integral representations of the solutions 


By means of the Laplace transformation and the Lorentz 
transformation we have reduced the time and the compressibility 
of the fluidity down to the two constants to* and x . The 
solution of our boundary value exercise can be represented 
generally by the following integrals over the lifting surface: 


G*(x.jr,z) = JJ w*(xV) K(x-x’ ,y-y* ,z) dx'dy' 
p*(x,y,z) = JJ\v*(x’,y') K-j ( x-x ' ,y-y ' , z ) dx'dy' 


( 10 ) 

(ID 


Because we have linearized our problem, the 
the differential nuclei of these integrals, 
can transform equation (11) into: 


nuclei K and are 
For this reason we 


P*( -x ' > - y ' > z ) - JJ w*(-x,-y) K.j (x-x ' ,y-y ' , z) ax dy . (12) 

From these equations one can easily derive the theory of the 
opposite flow formulated by A. H. Flax, as I have shown in 

another of my works. This theory runs in our case at the lifting 
surface 

z = 0 


JJ lw*(x,y) p*(x,y , + 0) - w*(x,y) p*(x,y, + 0)] dx dy = 0 . (13) 
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Here the horizontal line means opposite flow. Now we 
the following general statements for the down current 
and the pressure in direct and opposite flow: 


can make 

/41.8 


w *( x >y) = g( x »y) + f (x ,y ) ; p*( x ,y ) 
w *(*»y) * g(x.y) - f(x,y) ; P*(x,y) 


P g ( x -y) - P f (x,y) , 

P g (x,y) + Pf (x,y) . 


(14) 


Here g means even and f uneven down current dispersion 
the corresponding pressure dispersion. If „ e aubstitu 
(14) for (13), then we attain the tolerance condition 


t 


P g aIld P f 
e equation 


J7 t-(x,y ) p g (x,y) + g(x,y) p f (x,y)] dx dy - 0 


(15) 


If, for example, f(x,y) = x , g( x ,y) 
the tolerance condition of the lift 
equation (15) 


- 1 is set, then we obtain 
and torque coefficients from 


m + 
a 


0 , 


(16) 


to which I shall return later. All exact lifting surface 
theories must satisfy these tolerance conditions. 


The nuclei K and 
are, in spite of their 
complicated, as we shal 
hand, is the nucleus 
down current 


of the integral equations (10) and (11) 
different characteristic, quite 
1 see. Easily calculable, on the other 
of the inverse integral equation for the 
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"*(x,y,z) =2^// p*(x'.y',+o) M 1 ( x-x • ,y-y ' , z ) dx . dy . f 

( x >y 1 0) = /_“ explu*(a+x-Bv 4 2 +y 2 ) ] l+B^V^+y 2 ^ 

(oc 2 +y 2 )/a 2 +y 2 


( 17 ) 


(18) 


I a.erived the first approximation theory for oscillating lifting 
surfaces of finite wingspan from this integral equation in 1940 
[1. . According to this theory, one calculates the pressure 
distribution first on the assumption of even flow in all cross 
sections in accordance with the two-dimensional /41 . 9 

theory and then substitutes a new constant T (cj*,k,y) for the 
integral constant T («*). This new constant is calculated by 
solving a linear integral equation and contains the influence of 
the configuration of the lifting surface. The larger the span 
the smaller the difference of both functions. Later, E. Reissner 
derived this same theory by other means and it appears to be the 
sole three-dimensional unsteady lifting surface theory to date 
which has been employed to a wide extent in calculating 
directional oscillation flutter. 


In the steady boundary case cu * = 0 this theory extends into 
Prandt theory of the supporting line. None of these theories, 
however, satisfy the tolerance conditions (15) and (16), not even 
in the case of infinite span. The basic assumption of even flow 
is, therefore, inapplicable. Lift will be approximately 18* too 
high, and there also exist experimental indications for a mistake 
of this magnitude. The essential difference between a two- 
dimensional and a three-dimensional solution consists of the way 
a disturbance fades in infinity. From this vantage point it is 
surprising that the assumption of even flow in the case of three- 
dimensional problems does not create even greater mistakes. 
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The simple integral equati 
help us very much, because it 
surface which is given, but ra 
wish to calculate the pressure 
must invert it. This, however 
procedure. Generally such exe 
numerically and approximately, 
x and y with indefinite coeffi 
pressure at the lifting surfac 
(17), integrate and determine 
given down current is represen 


on (17), unfortunately, does not 
isn’t the pressure at the lifting 
ther at the down current. If we 
from this integral equation, we 
, is a difficult mathematical 
rcises can be solved only 

One could put up the polynomials 
cients for the desired /41 . 10 

e, introduce this into equation 
the coefficients such that the 
ted accurately. 


In the subsonic range, the desired counter pressure 
dispersion p* is uniform at the leading edge of the lifting 
surface, as I have mentioned already. Because only regular 
functions can be integrated numerically, we would have to know 
the manner of the singularity beforehand, separate from the start 
and integrate analytically. To date the same manner of 
singularity as in the case of even flow has been assumed without 
proof. 


In order to test the validity and the convergence of such 
numerical collocation methods, it is desirable to possess at 
least an exact solution of the three-dimensional unsteady lifting 
surface theory. Such a possibility exists, first off, only for 
the elliptic lifting surface, because this is the sole type of 
even lifting surface, for which separable solutions of the wave 
equation or of the Laplace equation exist in orthogonal curved 
coordinates. For this reason I want to present to you a theory 
of the oscillating elliptic lifting surface in the following. 
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5 . Theory of the oscillating elliptic lifting surface. 


5.1. Orthogonal coordinates and the Lame functions 

We shall consider an elliptic lifting surface with the 
designations used in Figure 1. The Cartesian coordinates x,y,z 
can be expressed by the ellipsoidal coordinates /41 . 1 1 

u,v,w as follows: 


x = cosh u cos v cos w 

y = k -1 \/ 1+k 2 sinh 2 u l/l -k 2 sin 2 v sin w 

r 2 p 

z = sinh u sin vV1-k' sin w 


/ 2 

Where k’ = v l-k . The square roots are always real and 
positive. In the boundary cases k = 0 and k’ = 0 these 
coordinates extend into those of the elliptic cylinder and of the 
spheroid. The surface u = 0 represents the elliptic lifting 
surface, upper side and under side. Its leading edge is given by 
u = v = 0, its trailing edge by u = 0, v = 7 T The Lame auxiliary 
quantities of these coordinates are 


U S; U <S; CD 
71 < V ^ 71 


S W £ 


(19) 


u = 


2 2 
1+k 9inh u 


o p ? 2 T" p 

(sinh u+sin v) ( 1+k sinh u-k' sin w) 

2 2 

1 -k sin^v 

(sinh u+sin v)(1-k sin 2 v-k' sin 2 vv) 

k 2 ( 1 -k* 2 sin 2 w) 

(1+k sinh u-k' sin w)(1-k sin v-k' 2 sin 2 w) 


( 20 ) 
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The surface element of the upper surface u = conal , is 


d o' 


dv dw 
V W 


( 21 ) 


Further we need the derivations according to z in the plane z = 0 


- U(0,v.,w) on the lifting surface ( u = o) , 

~5z = V(u,0,w) -A before the lifting surface (v = o) . 


( 22 ) 


We shall limit ourselves in the following to incompressible 

fluids ( H = 0) . For this reason we are seeking solutions of the 

Laplace equation A 2 <J> = 0, which are regular and periodic upon 

surfaces u — const, and disappear in infinity. These solutions 

consist of products of the Lame functions of the first and second 

order. The Lame f unctions of the first order are polynominal in 
. / 2 ! 2 ~ 

sm v, cos v and y 1-k sin v. For that reason they can /4l . 12 
be represented by associated Legendre spherical functions of the 
first order. For example, the Lame polynomial is 


Ec^(v) = sin v [l - Cj(k) sin 2 v] , 

“ Z P 3 ( cos v ) + ^ i5 - 4c^(k)] p|(cos v) , (23) 

Ecj(w) * l/l -k’ 2 sin 2 w [1 - c^(k') sin 2 v.-j 


The subscripts correspond to the nomenclature of Ince and 
Erdelyi; Ec“ are the functions which are even in cos v, Es m the 
uneven. The coefficients c™ are algebraic functions of k. The 
Lame functions of the second order are obtained from those of the 
first by insertion of the associated Legendre spherical functions 
of the second order (compare [2]): 
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FcJ(u) = C^(k)|{ Q^( cosh u) + ^ 15 - 4c^(k) ] Q^(cosh u)] ( 


(24) 


whenever u > 0 and k > 0 

There is still a second, more complicated representation of the 
Lamd functions of the second order using elliptic functions and 
integrals, which are also valid for u = 0 and k = 0. By using 
this representation one generally obtains the coefficients 


c >) - 1 + *>) K*(k) r 1 , (25) 

whereby a“ , b“ are algebraic functions and E’, K’ are complete 
elliptic normal integrations. 

I have normalized all Lamd functions, such that at the 
position zero either the function itself or its first derivation 
equals one. I designate the product from normalized Lame 
polynomials of the first order as ellipsoidal harmonics: 


Scj(v,w) = Ecj(v) Ec^(w) (26) 

A table of the first 24 Lame polynomials is to be found /41 . 13 
in the appendix of my printed report. 

In order to calculate with these functions, we must be able 
to diffentiate and integrate. On the lifting surface u = 0 the 
derivations of the ellipsoidal harmonics are 


Sx Sc 2n (v,v ') “ D(0,v,w) 
"5y Sc 2 n( v **) “ D(0,v,-w) 


I a° Ss 2r 


,.o . « • I 

l Scg - 1 


(27) 


1 
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Corresponding formulas are valid for the remaining ellipsoidal 
harmonics. For their integration the orthogonality relation is 
valid 


47i Vl+k 2 sinh^u 




(28) 

v.w) SCg(v.w) 

^ dv dw « 

r 0 - if 

m * r or 

n + s , 

1 . if 

m = r and 

n = s . 

the entire 

surface of 

the 


ellipsoid u = const. The coefficients 
functions of k. 




are algebraic 


5.2. Oscillating elliptic plates without approaching flow 

Now we shall consider the problem of the oscillating elliptic 
plate without approaching flow. By means of the Lame functions 
defined above we can make the following general statement for the 
regular velocity potential <P 


<p(u ,v,w) = k JJ <p u , (0,v' ,w' ) U( 0 , v ' ,w ' ) dor 1 • 

[ G ( u , v , w ; v ' , w ' ) + H(u,v f w;v' ,w' ) ] 


(29) 


Where = k/4n. For down current dispersions of the plate, which 
are even in y* or sin w* , we have the characteristic Green 


function 


/41. 14 
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CXD 

G ( u , v , w ; v ' ,w' ) - £ 
n = o 

n 

I 

m = o 

2x+1 

? c 2n + 1 

Fc 2m+1 (u) 
i 2n+1 v ; 

Sc 2 -!(v,vi) 

Sc 23+1 

bc 2n+1 

CD 

+ Z 

n 

Z 

YS 2m+2 

ys 2n+2 

r 2m+2 / \ 

Fs 2n+2^ u ^ 

S^v.w) 

Sa 2lE+2 

bs 2n+2 


n =. o m-0 


and for down current dispersions which are uneven in y* or sin w* 


H(u,v,w;v' ,w' ) * Z Z Y c 2n+2 fc 2n+2^ u ^ Sc 2n+2 (v,w ^ Sc 2n +2 (v/,v ^ 


2g+ 1 < 


2m+1 / f/ , 


naO m=o 


oo n 


(31) 


+ Z r Ta|s:| Fs^u) Sb-;5(v,w) (v». 


2n+2/ 


,._2a+2, 


n = o m = o 


At the plate u = 0 the given down current w* must be reproduced: 


V.’*( V ,w) 


Stp 


z=o 


9 u ( 0 ,v,w) U ( 0 , v , v; ) 


(32) 


Now the down current w* or the function upon the lifting 
surface is any one regular function from v and w and is, as such, 
representable in ellipsoidal harmonics using a progression. 
Because I normalized the functions Fc, Fs such that their 
derivations are 


Fc'(O) = Fb’(O) = 1 (33) 

and because the orthogonality relation equation (28) is valid, 
our statement (29) will be identically satisfied by the Green 
functions G and H in accord with equations (30) and (31). The 
first components of this statement for k’ = 0 are to be found in 
other representations already in the textbook by Lamb [3], Thus 
it is possible to call the statement (29) a classical statement. 

The derivations of the characteristic functions G and H 
according to z at the plate u = 0 are the so-called Dirac £- 
funct ions : 
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! 0, if v = v’ or w = w’ , 

(34) 

oo, if v = v* and w = w ’ , 

! 0, if v = v* or w = w’ , 

(35) 

oo, if v = v* or w = w* . 


5.3. Determination of the singular potential 

In order to attain a theory of the inflight lifting surfaces 
with Kutta conditions out of the classical statement (29), we 
must add a singular potential i p . Its derivation d*jj/dZ has to 
disappear upon the lifting surface u = 0, because we have already 
reproduced the down current w using the regular potential 
This requirement cannot, however, be completely satisfied, since 
periodic functions, in particular Fourier-series , do not permit 
representation of zero in the entire range of a period, in which 
all of the series coefficients do not disappear. If the singular 
potential is not to disappear identically, then its derivation 
must become infinite at least in one point of the lifting 
surface. 

We already know functions with this feature, which besides 
that satisfy the Laplace differential equation; their derivations 
are the k-functions seen in equations (34) and (35). If we wish 
to satisfy the Kutta condition of smooth air flow at the trailing 
edge, then we may lay the singular point at any one point v* = 0, 
w> = Wq of the leading edge. Then we obtain the singular 
velocity potential 

/41.16 
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( 36 ) 


♦ (u,v,w) =* \ Jf <p u ,(0,v< ,w') U(0,v’,w') da'- 

g( v ' , v* ' ) lid G( u, v ,\v; v ,w ) 0(0, v ,w ) 

V Q wO o o o o 

+ h(v',w') lid i!(u,v,w;v ,w ) U(0,v. ,w ) 

/ fl rO ^ u u 00 

Here g and h are functions which are still to be determined. In 
order to obtain flow at all points of the leading edge, that is, 
singular positions of the pressure, we can multiply equation (36) 
by any weight function 


a> co 

f(w 0 ) « Z a n cos n w 0 + r sin n w 0 (37) 

and integrate from w q = -c/2 to +c/2. Now we have 

Sc^(v,w) U(C,v,w) = 1 (38) 

Then the integrals of the individual series components become 


lie 


^ 2m+ 1 / \ 

jC 2r,+ 1 ( v ' w) 


X XUt I T 

V = 0 J-W/l 


lim ( 


Sc^(v,w) 

*/* sc|;+j(v,w) 


f(w) d W = A g:; = S c!-] tV ( 0 ,wi) , 

f(w) dw = b%:} - so^:i „ ..(o.wi) . 


(39) 


V=0 J-^ Sc](v fW ) ' "2n + 2 * -'*2n+2,v,w"'»"o / *■ (40) 


Because the functions g and h are still undetermined, we can fix 
the first coefficient of the series (37) arbitrarily: 


a Q = 1/c ; b 1 = 2/c 


(41) 


Then the coefficients become Aj = B 2 = 1. The expression to the 
right in the equations (39) and (40) is only a symbolic notation 
for the coefficients , B™ . Then in accordance with 
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equations (29), (36), (39) and (40) the entire velocity potential 
of the lifting surface becoies 

/41. 17 


<{> ( u , v , \v ; v I v ; 1 ) = <? + y = x(j /.•*( v',w' ) K(u,v,w;v',w' ) do ' , 
K ( u , v , w ; v ' .v ' ) = C(u,v f w;vjw' ) + g (v}w') G yI (u, v,.v;0 ,wj) 

+ HCu.v.srjvJw' ) + h ( v Jw ' ) H v , w ,(u,v,w;0 I wJ) 


(42) 


The constants contained in this statement can only be 
calculated later. First we must determine the functions g(v*,w’) 
and h(v’,w’). To this end we differentiate the nucleus K 
according to z on the plane z = 0 in front of the lifting 
surface, that is, for v = 0, and obtain 


as 

31 


= V(u,0,w) 


&K _ 1 

3v sir_h u 


lin K(u,v,ff;v!w') 
v=o 


0(0, v,w) 


(43) 


In order to reach the necessary symmetry of the nucleus in both 
the coordinates which are factored out and those which are not, 
we multiply equation (43) and integrate from w = -ff/2 to +tt/ 2. 
By solving linear differential equation (7) on the positive x- 
axis in front of the lifting surface, we obtain the additional 
condition 


GO 

K exp(-u* x)dx = 0 


(44) 


From the equations (42) to (44) there finally follows the nucleus 
of the velocity potential 
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( 45 ) 


K ( u , v , w ; v | w ' ) 
{ X 1 


+ X, 


\ exp(-u* coch u’) du' 

J ° i 

Gtu.v.wjv'w' ) , G^CuJO.w^jvJw' ) 

C y i (u,V,w;O t w^) . G vvl (u;0,^;0,wj) 
U(u,v,w;vjw‘ ) , 11 ^( 11 ', 0,-.v^;v;w' ) 

H v , w ,(u,v,w;0,wi) , H^^.UjO.w^O.v,*) 


) 


If K is to be a differential nucleus, then the following 
tolerance conditions with any one constant and Bust be 
valid: 


6 

b X 

Gfu.v.wjvJw' ) = 

C (0,0,' 

w i;v}w') , 

G , (u 

dx 

+ 3F-> 


0 

4 

v' 


G v (0,n, 

w 0 ;v;w’) , 

G v ,(u 

e 

37 

+ 3P° 

H(u,v,w;v;w' ) = 

H^(0,0 

. 0 ’ v ’ w ' ) * 

i 

H v’w' 



^(O.Tl 

• w 0 5 v ; w ')» 

H v'w* 


(46) 


Every differential nucleus satisfies the differential /41 . 18 

equation 


(|^ + ^|t) X(x-x’ ,y-y' ,z-z') = 0 


(47) 


From equations (7) and (47) -we obtain the nucleus of the integral 
representation of the pressure 
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-(u* + ^t)K 


( 48 ) 


K 


1 



We insert equation (45) into (48), integrate partially according 
to x and employ the tolerance condition (46). Then we finally 
obtain the integral representation of the pressure 


p*(u,v,w) = - \ jj w*(v;w’ ) K 1 (u,v,w;v;w' ) dcr ' , 

Kj (u,V,w; v}w 1 ) .= (u* + ^|r) G(u,v,w;v;w' ) + G^, (u,v,w; 0 ,'«*) • (49) 

•^-lG v (0,n l wJ;v;w , )T 1 - G T (OiO,w^T|w’)] 

+ (w* + ^p-) H(u,v ,w;v}w' ) + H Y , w ,(u t Y, W }0,wJ) • ( 59 ) 

• Cj-[H vw ( 0 , 7 t,w^;v;W)T 2 - H^CO.O.wJjv'w' )] • 


The constants T^ and Tg are given by 


T, = 


T~ = 


f ^ i i 

] ® vv i (u f O,w£}x,*r) exp(-u* oo3h u) du 

f 00 . j 

G vv , (u,0,w 0 ;0,w o ) exp(-u* cosh u)du 

•'O 

1 

fOO . . 

H v.vv'-.v ,(u ’ 0 ’V 0 ’ w o ) ex P(“ w * 003h u ) du 
Jo ' ’ J 


H ywv « w , (u,0,w^;ji,Wq) exp(-o)* cosh u) du 
>0 


(51) 


Now we must still determine the missing constants. We begin 
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According to this, the constants A and are independent 

of C 1 and Cg, while the constants A |a+| are proportional to C ]L 
and the constants B ^-2 are proportional to C 2 . If we insert 
these values into equation (46), then the constants Cj and C 2 
cancel out on both sides; for this reason they are not 
determinable by means of our statement. 

A further difficulty exists in the fact, that the set of 
equations (54) for n > 1 yields more conditions than are 

necessary for the calculation of the constants A For 

this reason all of the tolerance conditions cannot be satisfied. 

These difficulties probably can be avoided, if the extended 
statement for the nucleus of the velocity potential is made: 

/41.20 


G ( u , v , ; v ' , w 1 ) + H(u,v,w;v',w') 

r *r/z 

I c 

v'=o J -r/i- 


r< f.o l ^ C 2r + 1 Sc 2r+ 1 ^ v ' ' w ' ^ j r/i G ( u , v , w ; v', w*) U(0,v',w) ^^(W) dw 1 




(54a) 


By insertion of the derivative dK/bz for z = 0 in equation (44) 
and integration on skew lines y = const, one obtains for every 

Wichtung function f^ (w‘), f (w’), an integral 

equation in y 
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j, o exp(-a/x) dx Y(u,0,w) | Fc 2 °*](u) Sc^] t T ( 0,w) 

G v (u.O,w;v» U(0,v',wO f?®?](w') dw • } . 0 

v'ao-T/i v <;r+1 J 


Ho -5T/Z 

w »e arc sin ( . ■ - ) ; x * cosh u \J \ - 


kV 


l/l+k^sinh^u 


1+k sinh u 


-U k7 « i ; 


0 ; 


oo : u 


a> 


(54b) 


Solutions of these integral equations do not yet exist. It would 
be necessary to test if these solutions were unique and if the 
nucleus thusly calculated were a difference nucleus. 


The solution equation (50) obtained from the statement 
equation (42) is only approximately valid; it requires only the 
solution of finite linear equation sets, whose coefficients are 
easily calculable. The approximate solution equation (50) is 
applicable, although it is not unique. 


I am not of the opinion that the three-dimensional lifting 
surface theory in frictionless fluidity ought to be given up, 
because it isn’t unique. One could choose the constants and 
C 2 in such a way that measured and calculated aerodynamic 
coefficients agree as well as they can. Further, in boundary 
cases of infinite tension (k = 0), one could adapt them to the 
two-dimensional theory. The following constants correspond to 
this assumption: 


C, = - ycj Fc j (0) 
C 2 = - |C 2 Pc 2 (0) 


E’ Ik) 

’H' ~ 

(1 - 2^)2* + k X' 


(55) 
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With these constants I have calculated the aerodynamic 
coefficients, to which I shall return again. I cannot yet name 
any corresponding numerical values for the constants and Tg* 

On the basis of certain symmetric observations one can guess that 
the following boundary values are valid for stationary flow 
( CO* = 0) : 

if 


T 


1 



if 


k = 0 
k = 1 


(56) 


The solution obtained in equations (49) to (51) has the same 
form as those of the two-dimensional theory. It is valid for 
X = 0> can, however, be carried over in sections to 0, if the 
Lamd wave function is inserted in place of the Lam£ potential 
function. I carried out a corresponding transformation for the 
two-dimensional theory [4]. In the case of stationary flow ( GO * 

= 0) the solution obtained above is also valid for compressible 
fluids, because >( =tJ * B = 0, if one of the two factors 
disappears . 
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6 . Aerodynamic coefficients of the oscillating elliptic lifting 
surface. 


We now want to calculate an example showing the aerodynamic 
coefficients for an oscillating elliptic lifting surface with the 
bilinear down current dispersion 

w*(x,y) = a + bx + ck -1 y + ek -1 xy , 

w*(v,w) m [aSc^(v,w) + bSogt v ,w) + cSc^C v , w) ^ ^ 

+ eSsj(v,w)j U(0,v,w) 
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i 


i 

* 


i 

t 


The dimensionless aerodynamic coefficients, based on the elliptic 

2 

surface *1 /k and its half axis 1 and \/|< as lever arm, are 

for 


lift : 

pitching torque: 

rolling torque: 
deviation torque: 


k a,b 

_ k 

a 

11 p a,b (v ’ w ) 

Scj(v,w) 

U(0,v,w) 

d ff , 

m a,b 

k 

m. — 

71 

11 p a,b (v ’ w ) 

Ss 2 (v,w) 

U(0,v,w) 

dtr , 

n c,e 

k 

** — 

71 

11 p c,e (v ’ w ) 

Sc 2 (v,w) 

U ( 0 , v , w ) 

d<T , 

, e 

k 

« «— 
71 

W p c,e (v ' w ) 

Se|(v f w) 

U(0,v,w) 

dcr . 


(58) 


From equations (49), (50), (55), (57) and (58) we obtain the 
surprisingly simple pattern of the aerodynamic coefficients 


k Q = A[u>* +1 + T.J ] • n c ” +1+T 2 ] ’ 

~ k b “ “a “ A l‘?( 1+T i)] 5 -n e = q c = B [^(HT 2 )] , (59) 

E b = A («tu*+ 1(1-1^]; q e * B 100)*+ }(1-T 2 )] . 


This pattern is valid with other constants even for the two- 
dimensional theory and perhaps also for other non-elliptic 
lifting surface-contours. The constants in the two-dimensional 
theory are: 

/41.23 

A=|; ttA = ; B = 3 = 0 (60) 

and in the case of elliptic lifting surfaces 
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A 

cCA 

B 

SB 


F c ](0) 

* C 1 

^7 FaHco) 

Y3 2 

- -4 FcJ(O) 

yc 2 

P3^(0) 

ys, J 


4 1 

7 ^ 


17 

* 


1-k 

(2-k 2 )£' - k 2 K' 
1-k 2 

( 1-2k 2 )E' + k 2 K* 

,( 1^ 2 ) 2 


2( 1-k 2 +k 4 )E' - (k 2 +k 4 )K' 


1 


( 61 ) 


I would particularly like to point to the fact that the constants 
in equations (60) and (61) already arise from the classic theory 
of plates without oncoming flow oscillating in ideal fluidity. 
What is essentially new about the lifting surface theory is only 
the integration constants Tj and Tg, which can still be functions 
of the configuration of the lifting surface and of the reduced 
frequency co#, in the case of compressible fluidity as the 
Machnumber 0. 


7 • Numerical example for the boundary case k = 0 and k = 1 

In the boundary case of the elliptic lifting surface of 
infinite span (k - 0) and in the circular lifting surface (k - 1) 
there occur the constants contained in Table I 


Table I. Constants for k = 0 and k = 1. 


k 

E' 

K' 

A 

OA 

B 

BB 

*! 

A 2 

B 1 

$ 

y°i 

y® 2 

1 

y°2 

»•! 

0 

1 

1 

n 

7 

00 

n 

5 

4 

7 

8 

71 

2 

T7 

16 

T5 

$ 

2 

TCF 

64 ! 

1 

1 

1 

1 

1 

1 

1 

i 

3 

3 

15 

15 

15 

15 

105 

105 

4 bn 

1 57 5* 
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In order to compare available theories at least in the /41.24 

case of stationary flow, I have calculated the coefficients for k 
= 0 in accord with Prandtl’s hypothesis of even flow: 

k a = n ; k b = ‘ 7 5 

T. „ * . 

R c = 1 ; n e = “ B • 

For the circular lifting surface (k = 1) the numerical values 
from Krienes and Schade [5] are available. Further I have 
calculated the aerodynamic coefficients in accordance with 
equation (59) and Table I, whereby I inserted the values (56) for 
the T functions. The results are contained in Table II. 


n b = ° 


T5 


(62) 


Table II. Stationary aerodynamic coefficients 
for k = 0 and k = 1. 


k 

k a 

_' k b 

D a 

m b 

n c 

~ n e 


J q e 


0 

3,1416 

1,5708 

1,3333 

0 

0,7854 

0,3927 

0,2667 

■ 

| Procedure 
Prandt 1 


2,6667 

1,3333 

1,3333 

0 

0,5333 

0,2667 

0,2667 

1 

1 

Kiissner 

1 


J 

+ 17,8 

i 

+ 17,8 



0 

0 

+ 47,3 

+47,3 

0 

B 

1 

% Difference 











1 

0,8992 

» 

0,4718 

0,4659 

0,2191 

0,1276 

0,0582 

0,0554 

0,0299 

1 

i Krienes [5] 

: 


i) 

0,8488 

!i 

0,4244 

0,4244 

0,2122 

0,1132 

0,0566 

0,0566 

1 

0,0283 

| Kiissner 

, 

1 +i> ’ 9 | 

+ 11,2 

+ 9,8 1 

+3,3 

+ 12, 7| 

+2,8| 

-2,1 

*• i 

+6,0 

k 

% Difference 
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For the elliptic lifting surface of very large spans (k = 0) 
there are substantial differences between both theories at the 
lift and at rolling torque, while pitching torque and deviation 
torque agree. According to Prandtl the coefficients do not 
satisfy tolerance conditions 


k b + 



n + 
e 



(63) 


In the case of the circular lifting surface the differences are 
smaller; but the coefficients from Krienes do not satisfy the 
equation (63) either. The deviations represent 1.3X /41 . 25 

and 5. IX and are larger than the probable calculation errors. 


8 . Experimental testing and dispersion of the singularities. 

8.1. Displacement of the singularities at the edge of the 
lifting surface. 

It is not easy to test the lifting surface theory 
experimentally. Discrepancies occur as a result of the boundary 
layer and the finite airfoil thickness. If one uses very thin 
plates, discrepancies result from turbulence behind the leading 
edge. These deviations can be calculated only roughly or guessed 
at. The measurement values thusly corrected agree quite well 
with the new theory, less well with the old. 

I should like to mention a suggestion which was made by B.T. 
George, N. Rott and by me. In the case of oscillating wings in 
awater channel I have observed also that the trailing edge is 
weakly flowed around. The solution (50) with singularities on 
the leading edge could therefore be superimposed by a 
corresponding solution with singularities on the trailing edge. 
The components of both solutions are to be determined in such a 
manner that the measurement results are represented as well as 
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possible. If both portions are alike, then a solution with the 
stationary lift zero is obtained. 

A similar procedure is to be employed in the case of a 
thrusting lifting surface, whose flight direction is inclined at 
the angle a* against the symmetric plane y = 0. The Buler 
equation then becomes 

/41.26 

((J*- cos of ^ - Sin cc Jy) 4»*(x,y,z) + p*(x,y,z) = 0 . 

The leading edge of the lifting surface is shifted around the 
wing tips at the valve 

= arctg (k tg a *) (65) 

The singularities which are to be ordered along the new leading 
edge are once again to be determined such that the nuclei K and 
Kj become difference nuclei. I cannot discuss this any further 
here . 


8.2. Infinite pressure increase on the leading edge 

Much has been said about the infinite pressure increase near 
the leading edge of the lifting surface. Although the purpose of 
the lifting surface theory can only exist in the approximate 
calculation of the aerodynamic coefficients, this question has a 
certain interest. As long as no closed expression for the 
characteristic Green function is known, as in the two-dimensional 
case, the integral representation of the pressure equations (49) 
and (50) offer no answer to this question, because the series 
development of G for u = 0 becomes divergent. 
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One could try to calculate the pressure in a different manner 
out of the velocity potential. Corresponding to equation (48) 
one could, for example, differentiate the regular components of K 
according to (-x*), the singular ones according to x. One then 
obtains from equations (42) and (45) the pressure on the lifting 
surface u = 0 

P*(0,v,w' = -w *9 + f(v.w) + A U(0,v,w) [g 0 (v,w)+g 1 (v,*w)cos' v] (66) 

Here f and <p are finite amounts of ellipsoidal harmonics and A 
is a constant, all of which are dependent upon the given down 
current g Q and gj are even functions in cos v, which are 
independent of the down current. /41 . 27 

All of the statements of the three-dimensional lifting 
surface theory known to date are of the type given in equation 
(66) and this is the sole reason why I have presented this 
statement to you, because it is unacceptable. If we 
differentiate the function G v> (0,v,w;0,w^ ) using x, then we 
lose the characteristics, which are required. The Kutta 
condition is no longer satisfied and the singularity at the edges 
is no longer integratable. One can already be convinced of this 
fact in the two-dimensional lifting surface theory for which the 
Green function is 

00 i 

G(u,v,v') = - Y. n ex P ( -nu ) sin nv oin nv' , 

1 n (67) 

. 1 j n cosh u - cos ( v— v ' ) 

7 cosh u - cos(v+v') 
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Kinner and Krienes consider the functions g Q and g^, which are 
deduced from the Green function, to be freely selected with the 
side conditions, that p* of the Laplace equation suffices and 
that the Kutta condition and the down current condition be 
satisfied as best as possible. This is, however, only practical 
for a few points of the lifting surface. If the functions g Q 
(0,w) = gj (0,w) were final, then we would obtain the same kind 
of singularity of the pressure from equation (66) as in the case 
of even flow, that is, 1/sin v. Since, however, the statement 
(66) is unacceptable, this type of consideration has no 
convincing power. The oscillating plate without oncoming flow 
has a singularity like 1/sin v at the edge. 


9 . Solutions at any one temporal change of the down current. 

In closing I should like to discuss briefly the unsteady 
lifting surface theory at any one temporal change of the down 
current. At the outset we undertook a Laplace transformation of 
this general problem, in order to facilitate our task. After 
having solved the simplified task, we can cancel this 
transformation again by means of an inverse Laplace 
transformation. This occurs when we multiply the harmonic 
solution by the integral operator 


L 


1 fC + ieo 

2K±\c-ico- ex Pl“* 


(s-3 2 (s+x)) ) 



( 68 ) 
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If we employ this operator on the integral equations (17) and 
(18) for the down current, then we obtain the result 


J7~2 rs'=s 

vj(s,x,y) = — jff-JJdx'dy'l dp( s ,x,y ,+0) M 2 ( s-s ' ,x-x ' ,y-y ' ,0) , (69) 

Js'=o 


S + X+ 


■■ { s r +X) , v/enn s < 0 y/0s+x) 2 +y 2 /( 1-0 2 ) , 


M 0 (s t x,y ,0) =4 


y 2 ^s+x) 2 +y 2 /( 1-0 £ ) 


. x+fx 2 +y 2 


y 2 *£V 


, wenn s > s/( a+x) 2 +y 2 /( 1 -B 2 ) . 


(70) 


It is surprising, that the nucleus M 2 of this most common 
integral equation is such a simple algebraic function of the 
Cartesian coordinate. In order to eliminate the quadratic 
singularity of the nucleus M 2 , one must integrate equation (69) 
partially according to y’ and find the Cauchy principal value. 


A special, but equally interesting result is /41 . 29 

gained if we employ the integral operator (68) on the pressure of 
the harmonically oscillating elliptic lifting surface equations 
(49) and (50). Then we obtain the pressure in time point t = s 


p( 


u,v,w,s) = - \ JJ dcr^G+H)^! W(s,v',wO + w( 8t v»^|- (G+H) 

- G t (0,0,w£;v» w( 8 ,v»] (yjj 

+ H vV^ u,v,w;O ' w oHHy W (O,7i,w^;v / ,w0| 0 2 ( s-sOdif ( e\ v',*) 

* J s'=o 

- H^CO.O.w^v'w) w(8,v',w0] ]. 


The integrations constants T^ and T 2 are the picture functions of 
the object functions 
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^ rC - flco 

2xi J f _ ioo 




expu* a 


d u* 

u* 


if s > 0 (72) 


The function Uj(s) is to date only known in the two-dimesnional 
theory and was first calculated in 1924 by H. Wagner. 

We now consider the special wind dispersion 

w(e,x,y) = w(s+x,y) . (73) 

This down current occurs when an airplane flies into a stationary 

squall area. If we introduce this down current into equation 

(71), then only the singular components with the factors 1/C^ and 

l/Cg yield an amount different from zero to the pressure. For 

this reason there are in this case only two possible types of 

pressure dispersions at the lifting surface, a symmetrical 

proportional G , (0, v,w; 0 ,Wq) and an ant isymmetrical proportional 
V { 

H y , w , (0, v,w;0,w o ) . I determined the corresponding statement of 
the two-dimensional theory already in 1940. 
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The first twenty- four Lam6 polynomials 
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Ec °’ 2 

-c 0 ’ 2 

C > 


“ s 2 

Ea 1 * 3 

ta 3 

Es 1 ' 3 

tS 4 


Abbreviations : 


l£(v) 

i^) 

1 

d 

1 - c°’ 2 00 a 2 
dll - c°’ 2 (k) a 2 ] 

i 

a 

1 - c 2,0 (k' ) a 2 
a [ 1 - c|* 1 (k - ) a 2 3 

c 

c 

c d 

C 3 

■c [ 1 - a^dc) s 2 } 

c [ 1 - a 3,1 (k' ) a 2 J 

c d [ 1 - aj* 3 (k) a 2 ] 

c a[ 1 - aj* 2 (k’) a 2 ] 

S 

d 

s d 

d a 

a [ 1 - cj’ 3 (k) a 2 j 

d [ 1 - o 2,0 (k') a 2 3 

a d l 1 - cj’ 3 (k) a 2 ] 

dap- c 3,1 (k‘ ) a 2 j 

s c 

d c 

s c d 

d 2 3 

s c [ 1 - s| ,4 (k) a 2 3 

d c[1 - a 3,1 (k') a 2 ] 

a c d [ 1 - C 2 >2 (k) s 2 J 

d c a i 1 - c 2,0 (k’) a 2 ] 

• 


s = sin v 


sin w 


c = cos v 


cos w 


, ., 2.2 

d = 1-k sm 


l-k ,2 sin 2 w 
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c°* 2 (k) = i + k 2 ♦ i/T- k 2 +:-k 4 • 

c 3 ,2 (k) = 1 + 2k 2 ■+■ l/j" - k 2 4 4k 4 
s 3 ’ 5 0c) = 2 4 k 2 7 ^4 - k 2 4 k 4 
sj^(k) = 2 4 2k 2 7 /7 + k 2 + 4k 4 

C 3*>) * ? (2 + 2k 2 ? \f< - Ik? ♦ 4k 4 ) 

c]' 3 (k) - j (2 4 3k 2 ? y /7- 9k 2 4 9k 4 ) 

°4 ,4(k) e j (3 + 2 * 2 + ^9 - 9k 2 + 4k 4 ) 


i 


i 

l 
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